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We investigate in-plane quasiparticle tunneling across thin film grain boundary junctions (GBJs) 
of the electron-doped cuprate La2-i:Cea;Cu04, in magnetic fields up to B = 16 T, perpendicular 
to the Cu02 layers. The differential conductance in the superconducting state shows a zero bias 
conductance peak (ZBCP) due to zero energy surface Andreev bound states. With increasing 
temperature T, the ZBCP vanishes at the critical temperature « 29 K if B = 0, and at T = 12 K 
for B = 16 T. As the ZBCP is related to the macroscopic phase coherence of the superconducting 
state, we argue that the disappearance of the ZBCP at a field Bzbcp{T) must occur below the 
upper critical field Bc2{T) of the superconductor. We find Bzbcp{0) ~ 25T which is at least a 
factor of 2.5 higher than previous estimates of Bc2{0)- 
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Determining the magnetic field— temperature {B — T) 
phase diagram of high- Tc cuprates has been the focus of 
interest since the discovery of these materials. In con- 
trast to the case of conventional type II superconduc- 
tors, where the B — T phase diagram basically consists 
of the Meissner phase, the Shubnikov phase and the nor- 
mal state, the phase diagram of high- Tc cuprates is ex- 
tremely rich, exhibiting a variety of vortex phases and 
also a pseudogap region. However, particularly the tran- 
sition between the superconducting state and the normal 
state, and thus the relation between the superconducting 
and the pseudogap states, is hard to determine, not only 
due to the extremely large values of the upper critical 
field Bc2 in the case of hole-doped cuprates, but also be- 
cause of the presence of vortex liquid phases as well as 
strong fluctuation effects, leading to nonzero resistance in 
transport experiments long before Bc2 has been reached. 

In order to evaluate similarities and differences to 
their hole-doped counterparts, many investigations have 
been performed within the past years on electron 
doped materials like the single layer T' cuprates [ij 
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with composition L7i2-i^Cea;Cu04 {Ln = La, Pr, Nd). 
Pr2-:,Ce:,Cu04(PCCO) and Nd2-a;Ce^Cu04(NCCO) 
have been studied extensively, while there is less data 
published for La2-a;Ce2:Cu04(LCCO). The precise de- 
termination of the phase diagram of these materials is 
certainly of fundamental interest. 

In terms of Bc2, for PCCO and NCCO thin films near 
optimal doping, resistive measurements [2] or an analysis 
of the vortex Nernst signal [1, 0, @ revealed Bc2 (0) val- 
ues in the range of 7-10 T. For optimally doped LCCO 
an analysis of the vortex pinning strength yielded Bc2{0) 
around 9T However, the various methods applied 
to determine the upper critical field often yield incon- 
sistent results, see e. g. the discussion in In this 
paper we show that an analysis of Andreev bound states 
(ABS) causing a zero bias conductance peak (ZBCP) in 
the conductance spectra of cuprate grain boundary junc- 
tions (GBJs) yields a new lower bound for Bc2 which is 
at least a factor of 2.5 above previous estimates. 

ABS result from the constructive interference of An- 
dreev reflected electron like and hole like quasiparticles. 
If the quasiparticles experience a sign change of the su- 
perconducting order parameter upon reflection, the ABS 
appear at the Fermi energy, giving rise to the ZBCP 0. 
ABS-caused ZBCPs have been observed both in hole- 



2 



doped [i,i,[Tl [Hill and electron-doped [H, [Ti, [H, [3 

cuprates, where the sign change is due to the dx2-y2 sym- 
metry of the order parameter ITtI . This type of ZBCP 
can be observed with GBJs [llLlld| as well as for normal 
metal/superconductor junctions [1, 

or junctions between a conventional superconductor and 
a c?-wave superconductor . 

We note here that ZBCPs in the quasiparticle tunnel 
spectrum of superconductors can have other reasons than 
the formation of ABS, most prominent being the scat- 
tering of quasiparticles at magnetic impurities situated 
in the tunnel barrier, as described by the Appelbaum- 
Anderson model |20| . Such competing mechanisms 
can be identified, however, e. g. by analyzing the tem- 
perature and magnetic field depencence of the ZBCP. 

A ZBCP caused by ABS relies on the phase coherence 
of the elementary excitations above the Cooper pairing 
ground state. It thus should vanish when the macro- 
scopic phase coherence is lost, i. e., at the transition be- 
tween the superconducting state and the normal state. 
Such ZBCPs thus allow, at least in principle, to deter- 
mine Bc2 or at least to give a reasonable lower bound, 
provided that the ZBCPs do not vanish far below Bc2, 
e. g. due to the formation of a vortex liquid or other 
effects like a (field induced) shift of ABS to larger en- 
ergy. As we will see below, in the geometry used in our 
experiments an unsplitted ZBCP is clearly visible above 
the irreversibility line, allowing its use as a probe of the 
superconducting state in high magnetic fields. 

The 900 nm thin films for our study were deposited 
by molecular beam epitaxy with ozone as the oxidation 
gas dll on SrTiOa (STO) symmetric (ai = —a2, see 
inset of Fig. ^ [G01]-tilt bicrystal substrates with a to- 
tal misorientation angle (ai -I- a 2) of 24° and 30°. All 
films in this study had a Ce doping oi x ^ 0.08 (slightly 
underdoped, Tc ~ 29 K). The samples were patterned 
by photolithography and Ar ion millingwith junction 
widths ranging from 40 /im to 400 /im 'Slj. Transport 
data were collected with a current bias four point probe 
at 5 K< r < 40 K and magnetic fields B < 16T applied 
parallel to the thin film c-axis. d/(T^)/dl/-characteristics 
were measured with a lock-in amplifier (modulation fre- 
quency 1-5 kHz, lac = 1 — 10 /iA) and subsequently cal- 
ibrated with the numerical derivative of simultaneously 
recorded I{V)-cmves. Unless indicated otherwise, the 
data shown below are all for the same 24° GBJ. 

Fig. [1] shows resistivity p{T) of a bridge structure not 
crossing the grain boundary, for OT < B < lOT. Similar 
as for thin films of PCCO and LCCO [H, the on- 
set of Tc decreases with increasing B and the transition 
width is broadened. For i? > 6 T the resistance remains 
nonzero for T down to 5 K, and for larger fields we also 
observe an upturn in p{T), as already reported previously 
[2^ . [24j . We determine " transition temperatures" where 
the resistance of the (extrapolated) linear part of p(T) 
[cf. dashed line in Fig. [T] has dropped to, respectively 
10, 50 and 90 per cent of its normal state value, given 
by the intersection of the dashed line with the parabolic 
fit for p(T) above 35 K. The corresponding "transition 
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FIG. 1: (color online) p{T) of a LCCO thin film for OT < 
B < 10 T (in steps of IT). Dashed lines and full circles on 
the 6T curve indicate determination of temperatures where 
the resistance has dropped to, respectively 10, 50 and 90 per 
cent of its normal state value (grey line) . These values define 
the critical fields -Bio, S50 and _Bgo. The inset shows p(r) at 
B = up to 200 K with a parabolic fit (grey line). 



temperatures" or, respectively, transition fields -Bio, -B50 
and -Bgo , are shown further below in the B — T diagram 
of Fig. m Although the resistive transition broadens sig- 
nificantly, these fields do not differ strongly at a given 
temperature and tend to extrapolate to a zero tempera- 
ture value around 7T, suggesting that i?c2(0) is of this 
order. However, as we will discuss below, ZBCPs can be 
seen in much higher fields, suggesting e. g. superconduc- 
tivity below 12 K at 16 T. 

Fig. [2] (a) shows quasiparticle tunneling spectra (QTS) 
for a 24° GBJ, as obtained at 5K for B up to 16 T. At 
S = 0, a clear gap structure, that is a suppression in the 
quasiparticle density of states (QDOS), is observed, with 
symmetric coherence peaks at voltages Vgap — ±9mV. 
There is a small dip feature of unknown origin at volt- 
ages near ± 7 mV which disappears at higher fields. With 
increasing field the coherence peaks decrease and are 
strongly suppressed (although still present) for B > 6 T. 
For voltages well above Vgap the conductance increases 
linearly with voltage, cf. upper inset in (b), and is al- 
most independent of the applied field up to 6 T ^35| . For 
fields above 6 T the differential conductance strongly de- 
creases due to the onset of resistance in the film adding in 
series to the grain boundary resistance. This additional 
resistance has been subtracted in Fig. [2] (b), where one 
can see that the conductance in the subgap region in- 
creases monotonically with increasing field. The lower 
left inset in Fig. [2] (b) shows the integral over the tun- 
nel conductance, taken between -15 mV and 15 mV. The 
integral is field- independent, showing that the density of 
states is conserved. 

In the QTS of our samples we also observe a ZBCP, 
consistent with the previous study on optimally doped 
thin film LCCO GBJs [3. At T = 5K the ZBCP per- 
sists up to the highest fields achievable with our setup. 
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FIG. 2: (color online) Quasiparticle conductance spectra of 
a 24° LCCO thin film GB J at T = 5 K for T < B < 16 T 
in (steps of 2T). Graph (a) shows conductance as measured. 
In (b) the resistance of the leads appearing above 6 T has 
been subtracted. The inset in (a) shows an enlargement of 
the low bias region at 16 T. The upper inset in (b) shows the 
zero field conductance on a larger voltage scale at T = 5K 
(black line) and at T = 32 K (grey line), where the ZBCP has 
disappeared. The lower left inset shows the integral over the 
conductance curves taken for voltages between -15 mV and 
-1-15 mV. The lower right inset in (b) illustrates the geom- 
etry of the bicrystal GBJ. Vortices are indicated by circles, 
screening currents by arrows. 



i. e. 16 T. The ZBCP can already be seen clearly in the 
uncorrected data, as shown in the inset of Fig. [2] (a). We 
saw similar ZBCPs in 15 out of 19 samples both on 24° 
and 30° substrates. 

Fig. [3] shows the evolution of the ZBCP as a function 
of T at B = 16 T (a), and as a function of B at T = 13 K 
(b). Data are shown by thick black lines. The thin lines 
represent parabolic fits to the background quasiparticle 
conductance at low voltages. From the figures we see 
that, at i? = 16T the ZBCP disappears between 11 K 
and 12 K, while at T = 13K it disappears between 13 T 
and 14 T ^33^. 

As the ZBCP disappears reproducibly for a given T 
or B while increasing B or T, respectively, we can fol- 
low the disappearance of the ZBCP in a i? — T phase 
diagram, which is shown in Fig. [H Full symbols show 
the "critical" field Bzbcp where the ZBCP area Azbcp 
has reached zero for two different samples (squares: 30° 
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FIG. 3: (color online) Quasiparticle conductance dI{V)/dV 
of a 24° LCCO thin film GBJ in the vicinity of zero bias 
showing the disappearance of the ZBCP (a) at constant field 
B = 16 T for different T and (b) at constant temperature 
T = 13K for different B. Thin lines are parabolic fits to the 
quasiparticle conductance outside the ZBCP. 



GBJ, circles: 24° GBJ). We determined Azbcp by inte- 
grating the difference of the measured quasiparticle spec- 
tra and the parabolic background conductance [c. f. Fig. 
|3] in the range ± 3 mV. As AzBCp{B)\T=const and 
AzBCp{T)\B=const dccrcases nearly linearly with increas- 
ing B and T, respectively [s^], the intersection of the lin- 
ear fits with the Azbcp = axis defines Bzbcp (with 
vertical and horizontal error bars, respectively, in Fig. 
|4]). This procedure also allows to determine Bzbcp for 
B > IQT from AzBCp{B)\T=const (see inset of Fig. H]). 
Bzbcp increases monotonically with decreasing T, ex- 
trapolating to BzBCp{T = 0) « 25 T. The black and 
grey lines show Bc2{T), obtained from microscopic cal- 
culations (with Tc = 29K and Bc2(0) = 25 T), assuming 
either a 2D Fermi cylinder (black line) [IH] or a 3D Fermi 
sphere (grey line) [2^ . These lines may be considered as 
a lower bound for the true Bc2 of our samples. Note that 
Bc2{0) = 25T is at least a factor of 3.5 larger than the 
zero temperature extrapolation of Bgo- 

So far we have argued that the observation of the 
ZBCP extends the superconducting regime to magnetic 
fields that are much higher than anticipated. We next 
should discuss why we see a uniform ZBCP at all at large 
fields, without evidence for a splitting. A pronounced 
splitting would occur due to circulating currents if a sin- 
gle vortex were present near one side of the grain bound- 
271 . On the other hand, in the case of a symmetric 
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static vortex configuration, the screening currents that 
flow along the grain boundary are of equal strengths but 
opposite direction (cf. arrows in lower right inset in Fig. 
[^b). In this situation the shifts of the ABS (forming fast 
compared to the time scale of vortex fluctuations) should 
cancel and no splitting is expected. In the case of a vor- 
tex liquid, which is likely to be present in the high field 
region we discuss here, the screening currents on the two 
sides of the GBJ will fluctuate, but are likely to have zero 
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FIG. 4: (color online) Magnetic field vs. temperature phase 
diagram showing the fields Bio (open squares), B50 (open 
diamonds) and Ego (open triangles) , as determined from p{T) 
curves of Fig. [T] together with the field Bzbcp where the 
ZBCP in the quasiparticle tunnel spectrum disappears. Full 
circles: 24° GBJ, fuU squares: 30° GBJ. The black and the 
grey lines correspond to Bc2{T), calculated for a 2D Fermi 
cylinder and a 3D Fermi sphere, respectively. The inset shows 
the ZBCP area vs. B for different temperatures. Lines are 
linear fits, which, by extrapolation to zero ZBCP area, define 
values of Bzbcp for B > 16T. 

average. Thus again the energy shifts of the ABS will be 
washed out during the time of measurement, naturally 
explaining the experimental observation. On the other 



hand, at lower fields, where static but disordered vortex 
lattices may form, some net screening current is flowing 
along the GBJ. In this case the ABS at zero energy should 
undergo a Doppler splitting [27^ . Indeed we observe that 
the ZBCP becomes structured and often asymmetric in 
fields below the irreversibility line, possibly allowing to 
use the ZBCP to probe details of the vortex lattices. A 
detailed discussion, however, is beyond the scope of this 
paper and will be presented somewhere else [28| . 

Returning to the QTS as shown in Fig. [21 we note that 
the coherence peaks - often used as an indicator for the 
superconducting state - are suppressed already at 50 per 
cent of Bc2 or so. This observation is indeed consistent 
with calculations based on the quasiclassical Eilenberger 
equations of the density of states in the vortex state, 
as given in Fig. 5 of Ref. where, at 0.5i3c2 the 

coherence peaks are hardly visible. 

In conclusion, we have shown that a zero bias con- 
ductance peak can be seen in the quasiparticle tunneling 
spectra of LCCO grain boundary junctions at magnetic 
fields that are far above previous estimates of Bc2- Par- 
ticularly, a zero temperature extrapolation suggests that 
Sc2(0) is at least 25 T. Extending superconductivity to 
such high fields shrinks the region where a pseudogap 
phase may exist. 
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